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TEARING  INSTABILITY  IN  AN  ANISOTROPIC  NEUTRAL  SHEET 

I.  INTRODUCTION 

In  this  paper,  we  investigate  the  tearing-mode  stability  properties  of 

a  field-reversed  plasma  layer  whose  temperature  distribution  is  anisotropic 

(T  *  T^).  T^e  symbols  1  and  ^  refer  to  directions  perpendicular  and 

parallel  to  the  equilibrium  magnetic  field,  respectively.  Figure  1  shows  a 

schematic  drawing  of  the  geometry  and  the  coordinate  system.  The  equili- 

o  A 

brium  magnetic  field  B^z)  =  indicated  by  the  solid  lines,  is 

generated  self-consistently  by  the  current  J^(z)  =  J0(z)y^  with  no  external 
magnetic  field.  The  magnetic  field  reverses  its  direction  in  the  plane 
z  -  0  and  JQ(z)  is  symmetric  in  z.  Both  ions  and  electrons  are  assumed  to 
contribute  to  the  equilibrium  current.  The  number  density  nQ(z)  is  also 
shown  in  the  figure.  In  this  one-dimensional  equilibrium,  all  physical 

quantities  depend  on  z  alone.  The  plasma  is  assumed  to  consist  of 

collisionless  ions  and  electrons  so  that  Vlasov  equations  are  used  for  both 
species. 

It  is  well-known  that  inertia  of  the  current  carriers  can  lead  to 

collisionless  tearing  instability  in  a  system  such  as  the  one  described 

above.  Considerable  work  has  been  done  on  the  basic  collisionless  tearing- 

1—8 

mode  properties  of  the  neutral  sheet  configuration  .  An  example  of 

physical  systems  that  may  be  modelled  by  this  configuration  is  the  earth's 
9 

magnetotail  .  Subsequent  to  the  early  work,  tearing-mode  results  have  been 
considered  in  connection  with  the  magnetotail10,11.  In  these  studies,  the 
neutral  sheets  are  generally  assumed  to  have  isotropic  temperature 
distributions  (T^  =  T^).  However,  in  a  collisionless  plasma,  the  motion  of 
particles  parallel  to  the  magnetic  field  is  decoupled  from  the 

perpendicular  motion  and  temperature  anisotropy  can  be  maintained  even  in 

12 

thermal  equilibrium.  Laval  and  Pellat  used  an  energy  principle  analysis 
to  show  that  collisionless  tearing-mode  properties  can  be  strongly  modified 
by  weak  temperature  anisotropy.  In  particular,  it  was  found  that  the 
k  |  mode  can  be  stabilized  by  a  very  small  degree  of  electron 
temperature  anisotropy:  T  /T  <  (1-a  /5)  where  a_  is  the  electron 

“1  “P  6  “ 

gyroradius  and  6  is  the  layer  half-thickness.  In  this  work,  neither  the 
dispersion  relation  nor  the  marginal  condition  for  the  anisotropic  case  was 
given.  Recently,  Chen  and  Davidson13  carried  out  a  Vlasov-fluid  analysis 
of  a  field-rfversed  ion  layer  at  marginal  stability  using  approximate 
Manuacript  approved  July  14,  1983. 
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orbits  in  a  cylindrical  coordinate  system.  It  was  found  that  relatively 
small  temperature  anisotropy  in  ions  significantly  modified  the  stability 
boundaries.  Consistent  with  the  results  of  Ref.  12,  the  range  of  unstable 

wavenumbers  is  reduced  for  T.  ,/T.  <  1  and  is  increased  for  T  /T .  >  1. 

More  recently,  Basu  and  Coppi1  studied  a  fluid-like  "field  swelling" 

instability  in  an  anisotropic  plasma  in  which  the  magnetic  pressure  is 
comparable  to  the  particle  pressure  and  in  which  there  is  no  equilibrium 
current.  This  analysis  is  local  and  is  based  on  the  fluid  equations. 

However,  in  the  neutral  sheet  configuration,  the  instability  is 
intrinsically  nonlocal  and  the  plasma  0  is  much  greater  than  unity  in  the 
field-null  region  so  that  the  fluid  equations  are  invalid.  Thus,  the 
treatment  of  Ref.  14  is  not  applicable  to  the  tearing-mode  instability 
which  is  due  to  resonant  particles  crossing  the  field-null  region. 

A  difficulty  in  treating  the  neutral  sheet  configuration  is  the 
complicated  particle  orbits*  In  order  to  make  the  problem  tractable 
analytically,  "straight-line"  approximations  have  often  been  used  ’  for 
the  particle  orbits  crossing  the  null  plane.  Numerically,  the  problem  is 
also  difficult.  In  a  noteworthy  but  rarely  referenced  work,  Holdren18 
utilized  an  integrodif f erential  equation  formalism  to  show  the  feasibility 
to  treat  the  various  orbits  exactly.  In  order  to  minimize  the  numerical 
errors,  relativistic  electrons  with  large  orbits  were  used  and  no 
definitive  comparison  can  be  made  between  the  relativistic  numerical 
results  and  the  previous  nonrelativistic,  approximate  results'*.  In 
addition,  the  numerics  required  are  prohibitive  and  it  is  in  general 
difficult  to  use  this  method. 

In  the  present  model,  we  will  adopt  the  straight-line  orbits  but 
improve  the  treatment  of  the  ion  orbits.  The  effects  of  the  small  Larmor 
radius  orbits  will  also  be  considered.  The  anisotropic  collisionless 
tearing-mode  analysis  is  carried  out  for  perturbations  of  the  form 

A 

4>(x, z , t )  =  <|i(z)exp(ikx-iut)  where  the  wavevector  k_  is  taken  to  be  parallel 
to  the  equilibrium  magnetic  field.  In  Fig.  1,  the  dashed  lines  indicate 
schematically  the  magnetic  field  lines  including  the  mirror-like 
perturbation  given  above.  As  a  general  remark,  the  problem  treated  here  is 
different  from  that  of  tearing  instabilities  in  tokamak  discharges18 
because  the  tokamak  geometry  would  have  an  applied  magnetic  field  in  the  y- 
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direction  (Fig.  1).  This  alters  the  prticle  orbits  significantly.  The 
configuration  studied  here  is  more  closely  related  to  certain  0-pinches  and 
ion  layers. 

In  Section  II,  we  present  the  basic  linear  analysis  of  the 
perturbation.  In  Section  III,  simplified  particle  orbits  are  used  to 
obtain  the  eigenvalue  equation.  This  equation  is  solved  to  give  the  linear 
anisotropic  dispersion  relation.  Section  IV  discusses  the  results. 

II.  FORMULATION 

Figure  1  shows  a  schematic  drawing  of  a  neutral  sheet.  We  assume  that 
there  is  no  equilibrium  electrostatic  field  and  no  bulk  motion  of  the 
plasma  (i.e.,  E_  x  B_  ■  0).  We  consider  a  class  of  equilibria  described  by 
foj  *  foj<<Hij  '  V jpy j »  H„j)  where  J  "  e,i  and  (H^,  Pyj,  H^)  are  the 
single-particle  constants  of  the  motion  in  the  equilibrium  magnetic 
field.  Here,  -  (nij/2)(v2  +  v2),  Pyj  -  mjvy  +  (qj/c)A°(z)  , 

H  “  (m  /2)v2  and  A°(z)  Is  the  vector  potential  for  the  equilibrium 

"  j  j  ^  y 

magnetic  field,  where  qj  is  the  electric  charge,  mj  is  the  mass,  and  = 
constant  is  the  mean  drift  velocity  of  the  species  j. 

A.  Linear  Eigenvalue  Equation 

The  analysis  is  carried  out  for  perturbations  of  the  form  <|;(x,z,t)  = 

A 

<|»(z)exp(ikx-iuit)  where  the  wavevector  k^  -  W  is  parallel  to  the  equilibrium 
magnetic  field.  The  magnetic  field  configuration  including  this 
perturbation  is  indicated  by  dashed  lines  in  Fig.  1.  The  characteristic 
frequency  is  low  with  lu>l  <<  u>cl  where  (ocl  -  eBQ/m1c  is  the  ion  cyclotron 
frequency  in  the  asymptotic  magnetic  field  Bq  =  B°(z  ■  »)•  Moreover,  we 
consider  the  case  where  kvg  ~  u>c^  which  is  true  for  a  wide  range  of 
parameters  of  interest  so  that  we  also  have  typically  | col  «  kvg,  where  vg 
is  the  electron  thermal  velocity.  This  will  be  verified  a  posteriori. 

Using  the  standard  method  of  characteristics,  the  first-order  Vlasov 
distribution  function  for  each  species  is 


fi  4  *  "  -  I  dt'(E.  +  -  v'  x  B.  )  .  f  , 
lj  m  1  —1  c  —  —1'  3v'  o' 
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where  (it',  v")  represents  the  particle  orbits  in  the  equilibrium  field  with 
the  conditions  x'(t'  m  t)  ■  x»  *  O  “  v.  anc*  ♦(t  *  ““)  ”  0  for  all 
perturbation  quantities.  The  time  integration  is  carried  out  along  the 
unperturbed  orbits.  In  the  above  equation  and  in  the  remainder  of  the 
paper,  the  species  index  j  «  e,i  is  suppressed  where  no  confusion  arises. 
The  perturbed  fields  are  given  by 


”  1  X  A_1  , 

(2) 

1 

^’cF  * 

(3) 

B  *  ^  J 

S-i  ~-i  * 

(4) 

and 


V  .  E t  -  4irPl, 


(5) 


where  and  $  are  the  perturbed  vector  and  scalar  potentials.  In  this 

paper,  we  choose  the  Coulomb  gauge  *  0).  Because  the 

characteristic  frequency  of  the  perturbation  is  low,  w  <<  to  . ,  we  assume 

Cl 

charge  neutrality  to  first  order.  Note  that  Bly(x,z,t)  *  0  for  the  mirror- 
like  perturbations  so  that  »  Jly(x,z,t)£_ and  A^  *  i)i(x,z,t)y^ 


For  the  general  form  fQ 

3f_  -  3f 

3v  W 


fo  ~  vjpy.  V>  n°te  that 


_  -  3f  3f 

— —  *  m( v  —  Vy)  — —  +  mf -  —  — — v  x. 

av  “'■X  L >  ail  “'-3H  3H,'  x- 
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Using  Eqs.  (2)  and  (3)  in  Eq.  (1),  straightforward  calculation  yields 

3f  g  f  gf 

fu  ‘  "i  mf  '  Ba*>  +  luS,]  •  lk<1iC3H7  -  3ir>Vy 

1  -LI 


(6) 


where  Sy  is  the  orbit  integral  along  the  unperturbed  orbits  given  by 


S 

y 


/  dt'($ 


-  —  V">) 

c  yY 


(7) 


and 
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a 


is  the  uniform  mean  drift  velocity  normalized  to  the  speed  of  light  c 
Here,  use  has  been  made  of  3/3t  +  -ini  and  3/3x  +  ik.  For  the  isotropic 
case,  the  last  term  in  Eq.  (6)  vanishes  and  we  trivially  recover  the 
previous  results®.  For  the  anisotropic  case,  estimating  vx  by  the  typical 
thermal  velocity  vth,  we  note  that  the  ratio  of  the  second  term  to  the 
third  term  is  of  the  order  of  w/(kvth)  which  makes  the  last  term  in  Eq.  (6) 
small  in  comparison  with  the  first  two  terms  for  low  frequency 
perturbations  except  for  a  very  small  degree  of  anisotropy.  This  means 
that  the  marginal  state13  and  the  growth  rate  can  be  significantly  modified 
in  comparison  with  the  isotropic  case.  It  is  the  purpose  of  this  paper  to 
give  a  quantitative  illustration  of  this  point. 

Another  consequence  of  the  low-frequency  nature  of  the  perturbation  is 

g 

that  the  perturbed  scalar  potential  <f>  is  much  smaller  than  (vy/c)i|>  .  This 

can  be  seen  by  noting  that  the  electric  field  Elx  along  the  magnetic  field 
is  at  most  comparable  to  the  components  perpendicular  to  the  magnetic 
field,  say,  Ely.  That  is  to  say,  Elx  <  Ely.  This  implies 


For  the  low  frequency  case  (u>  <<  kv^),  Is  ranch  less  than  (vy/c)tj;  where 
vy  -  vth  for  the  typical  thermal  particle.  Physically,  this  is  simply  a 
statement  of  quasi-neutrality  at  low  frequencies.  Although  it  is 
theoretically  straightforward  to  include  <j>,  we  will  adopt  the 
simplification  of  neglecting  $  in  comparison  with  (vy/c)i|i  .  Then,  the 
approximate  perturbed  distribution  function  is 


f 


13 


3f  3f0  3fQ 

•q®  wL  ♦ +  i0)q  w~L  Vikqf  ■ 


)v  S 

*  V  ' 


x  y 


(8) 


where  the  orbit  integral  [Eq.  (7)]  has  been  simplified  to 


S 


y 


j. 

c 


t 


(9) 


The  firr  term  in  'Jq.  (8)  is  the  usual  adiabatic  term  which  does  not  depend 
on  the  le.  particle  orbits.  The  second  term  represents  the  change  in 
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momentum  and  the  third  term  represents  the  change  in  energy  of  the 
particles.  In  particular,  ikqSy  is  the  time-integrated  Lorentz  force 
acting  on  the  particles  in  the  x-direction. 

B.  Anisotropic  Plasma  Layer  in  Thermal  Equilibrium 


In  the  preceding  section,  we  have  described  the  general  theoretical 
framework  appropriate  for  equilibria  of  the  type  fQ ^ (H^  -  V^P^,H^).  To 
examine  the  stability  properties  in  detail,  we  specialize  to  the  case  where 
the  plasma  layer  consists  of  two  species  (electrons  and  ions),  both 
satisfying  the  two-temperature  Maxwellian  distribution  given  by 


=  ^7^7 exp(-  £ ‘  ’77 Hu>-  <10) 

where  =  (m.j/2)(v2  +  V2),  =  mjVj  +  (q j/c)A°(z) ,  =  (m./2)V2, 

V.  =  constant  is  the  mean  velocity  of  the  species  j,  and  A°(z)  is  the 

^  y 

vector  potential  for  the  equilibrium  magnetic  field.  As  stated  before,  the 
electrostatic  field  is  taken  to  be  zero  in  the  frame  of  the  layer. 


The  equilibrium  particle  density  and  magnetic  field  profiles 
corresponding  to  Eq.  (10)  depend  only  on  the  perpendicular  temperature  and 
are  well-known1 ^ 


n  (z)  =  n  sech2  (|) , 


(H) 


and 


2cTx/6 


Vz)  ~  e(V~=V  )  tanh  (j)' 
i  e 


(12) 


where  6  is  the  characteristic  half-thickness  of  the  plasma  layer  given  by 

c2T 


,  -  r  i)1/2  1 

6  =  S  A  ^  e(V  -V  ) 
2irn  i  e 


(13) 


Here,  T^  =  T^^  +  TfiX>  and  nQ  is  the  maximum  number  density  at  z  *  0  for 


both  species.  For  z  +  <»,  we  have  B  (z  +  •)  ■  B  with 

x  o 


B  =  /87m  T  , 
0  o  1 


(U) 
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where  8Q  >  0  is  chosen  without  loss  of  generality.  This  choice  implies 
(Vi-V.)  >  0  so  that  the  total  current  is  in  the  +y-direction,  given  by 

J^z)  =*  e(V1-Ve)nQ(z)y.  (15) 

Because  of  charge  neutrality,  the  equilibrium  ion  and  electron  densities 
are  equal.  It  is  easy  to  show  that  the  zero  electrostatic  field  condition 
is  equivalent  to 


T  ,V 
il  e 


+  T  V. 
el  i 


0. 


(16) 


If 


we  denote  the  gyroradius  of  the  thermal  particles  by 


a j  s  ”her<!  vn  5  /2V*j 


and  u) 


eB  /m.c  with  B„  given  by 
o  J  °  6  y 


cj  "V"T 

Eq.  (14),  then  it  is  easy  to  show  by  using  Eqs.  (13)  and  (16)  that 


6  vjl 


(17) 


and 


Vi 

CTjl 


V 


(18) 


As  a  general  remark,  the  equilibrium  configuration  described  above  is 
obtained  by  balancing  the  Lorentz  force  c  *(J^  x  B^)  with  the  perpendicular 
pressure  gradient  ~(3/3z) [nQ(z)I  ] . 


III.  Stability  Properties  for  an  Anisotropic  Neutral  Sheet 

A.  Model  Orbits 

In  order  to  determine  the  stability  properties  of  the  system  described 
by  Eqs.  (10)-(18),  the  orbit  integral  Sy  [Eq.  (9)]  must  be  evaluated. 
Generally,  the  orbits  in  the  equilibrium  field  are  cotnplicted. ^  In  the 
previous  calculations  on  the  isotropic  neutral  sheets,  "straight-line" 
orbits  ’  have  been  used  to  evaluate  the  orbit  integral.  In  these 
approximations,  the  neutral  sheet  is  divided  into  two  regions.  In  the 
inner  region,  |z|  <  ~  /a^d  where  the  magnetic  field  is  weak,  the  orbits 
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are  taken  to  be  straight  lines  across  the  null-plane  (z  »  0)  and  reflected 


from  the  surfaces  at  |z|  a  d 


Thus,  Vy  is  nearly  constant  in  this 


region.  In  addition,  the  perturbed  vector  potentrial  -  A^  is  assumed 
to  be  constant  in  the  inner  region  (constant-^  approximation).  In  the 
outer  region,  |zl  >  d^,  the  particles  are  assumed  to  execute  nonaxis¬ 
crossing  small  Larraor  radius  orbits  and  these  orbits  are  neglected.  In 
particular,  in  Ref.  8,  the  dispersion  relation  is  obtained  by  matching  the 
inner  and  outer  solutions  at  z  =  /2ae<5.  This  approximation  neglects  the 
axis-crossing  ion  orbits  that  extend  far  beyond  the  electron  inner  region 
since  d^/dg  ~  (m^/nig)1^.  Although  the  effects  of  the  axis-crossing 
electrons  are  greater  than  those  of  the  axis-crossing  ions,  the  ion  effects 
dominate  in  the  intermediate  region  dg  <  |zl  <  d^,  as  will  be  shown  in  this 
section.  In  the  isotropic  case,  it  will  turn  out  that  the  neglect  of  the 
intermediate  region  does  not  change  the  dispersion  relation 
substantially.  However,  in  the  anisotropic  case,  the  ion  orbits  in  the 
intermediate  region  affects  the  dispersion  relation  significantly.  In  the 
present  analysis,  we  include  the  three  regions  (Fig.  2)  and  define  the 
boundary  surfaces  at 


d  =  /2a  5 
e  e 


(19) 


for  the  electrons  and 


dj  = 


(20) 


for  the  ions,  where  Sj  *  vii^u)c1  ls  tlle  ^arinor  radfus  of  a  thermal  particle 
with  »  /2T •  Physically,  dg  is  the  distance  where  the  local 
electron  Larmor  radius  ag(dg)  =  ve^/u»cg(de)  Is  equal  to  dg/2. 

Here,  <*>ce(de)  “  eB°(de)/mec.  d^  is  the  distance  where  the  local  ion 
Larmor  radius  aj(d^)  =  v^/u>ci(d^)  is  equal  to  d^ .  The  factor  of  /T  is 
somewhat  arbitrary  and  is  chosen  to  avoid  overestimating  the  large  ion 
orbits.  Moreover,  this  choice  allows  the  use  of  the  constant¬ 

ly  approximation  in  the  region  |z|/5  <  d^S  <<  1. 


The  orbit  integral  Sy  [Eq.  (9)]  will  now  be  evaluated  for  the  three 
regions  as  shown  in  Fig.  2.  In  the  respective  inner  region  for  each 
species,  the  orbits  consists  of  nearly  straight  segments8’10  and  vy  is 
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In  deriving  Eqs.  (22)  and  (23),  use  has  been  made  of  the  equilibrium 
properties,  Eqs.  (13),  (14),  (16)  and  (18),  and  ■  /2T 

In  the  outer  regions,  regions  II  and  III  for  electrons  and  region  III 
for  ions,  the  particles  mainly  execute  VB  drift  motion  with  the  drift 


11 


velocity  VjJ  given  by 


Jj 


a  ®cv2 

"^V 


slnh  2(|), 


(25) 


where  v2  =  v2  +  v2  is  a  constant  of  the  motion.  This  expression  is  a  good 
1  '  y  z 

approximation  for  the  guiding  center  motion  in  the  outer  regions  for  most 
particles  (a^  «  5)  except  for  the  ones  with  kinetic  energy  much  greater 

than  the  thermal  energy.  Note  that,  for  thermal  particles  (v  =  v,  ), 

*  —2  1  ^  ^ 

Vp/v^  is  of  the  order  of  (a.j/6)  sinh  (z/6)  so  that  VjJ/v^  is  of  order 

unity  near  z  =  dy  The  actual  motion  of  a  particle  is  then 


Vp  +  v  cos(iocjt  + 


X) 


and 


v  »  v  sin(a)  ,t  +  X), 
z  cj 

where  w  .  is  the  local  cyclotron  frequency,  v  is  the  gyration  velocity  in 
^  J 

the  VD  frame  and  X  is  the  phase  angle.  Then,  the  approximate  orbit 
integral  Sy  for  the  outer  regions  is 

Sout  *  ic  lvp(kvx~<‘>)  1  v|»  ~  c  1  v  /  dt'cos(u>c .  t"  +  X)i|». 
y  — oo  ^ 

The  second  term  represents  the  oscillatory  motion  in  the  guiding  center 
frame.  Taking  the  phase  average  over  X  for  the  low  frequency  perturbation, 
it  is  easy  to  see  that  the  second  term  averages  to  zero  and  we  have 


S°ut  -  ic  1V^(kvx“o))  1\|»(x,z,t)  . 

Integrating  over  the  velocity  space  and  after  some  algebra,  we  obtain 


(26) 


out 

li 


c 


4im  62 


t1  -  (dH 


ii 


z(^kvtl) 

sinh2(z/$) 


(27) 


W(w/kv  ) 

+  a.  -  ]n  (z)t|/, 

sinh2(z/S) 


and 


.out 

le 


2T  .  Z(w/kv  .) 

c _ el  r ,  _  /  oi  \  _ el 

*  „  T  l 1  (w  7 

4nn  62  1  el  sinh2(z/6) 


+  o_ 


W(u/kve)|) 

sinh2(z/6) 


]n  <zH- 


(28) 


In  obtaining  these  expressions,  use  has  been  made  of  the  equilibrium 
relationships  Eqs.  (17)  and  (18),  and  a2/62  «  1  has  been  dropped. 

The  third  term  in  the  square  brackets  of  each  equation  is  the 
anisotropy  term.  In  Ref.  8,  where  the  isotropic  case  is  treated,  the 
second  term  in  each  equation  is  neglected  by  using  the  ordering 

(w/kv  )  -  (a,/6)^^  <<  1  and  the  ion  contribution  J*!*  is  neglected  in 
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region  II.  Galeev  and  Zelenyi  estimated  the  drift  velocity  by  [Eq. 

(25))  without  the  factor  sinh~^(z/6)  and  concluded  that  this  contribution 

is  small.  However,  as  the  discussion  following  Eq.  (25)  shows,  this  term 

can  be  of  order  unity  in  Eqs.  (27)  and  (38)  near  z  ■  dg.  Moreover, 

Q  1  A 

Dobrowolny  and  Galeev  and  Zelenyi1  matched  the  interior  (|z|  <  dg)  and 

exterior  (|z|  >  dj)  solutions  at  z  ■  dg  by  arguing  that  the  interior 

electron  contribution  to  the  perturbed  current  density  [Eq.  (23)]  is  much 

greater  than  the  external  ion  contribution  [Eq.  (27)]  for  Izl  >  dg.  This 

argument  overlooks  the  fact  that  the  axis-crossing  ion  orbits  extend  far 

beyond  the  electron  inner  region  (|z|  -  dg)  so  that  the  relevant  comparison 

in  the  intermediate  region  (region  II)  is  between  the  exterior  electrons 

[Eq.  (28)]  and  interior  ions  [Eq.  (22)].  For  the  case  in  which  ions  and 

electrons  are  both  anisotropic,  the  ratio  at  z  ■  dg  in  region  II 

is  approximately  (a^/ae)(me/m^) (6/ae)  which  is  typicality  of  order 

unity.  For  the  case  in  which  only  ions  are  anisotropic  (a  *  0),  the 

5/2  e 

ratio  is  approximately  a, (6/ae)  (me/tn^)  which  is  much  greater  than 

unity.  As  z  increases  to  d^ ,  the  ratio  increases  because  of  the  factor 
-2  — 

sinh  (z)  in  the  electron  contribution.  Thus,  in  general,  the  ion 
contribution  is  not  negligible  in  the  intermediate  region.  Moreover,  this 
Intermediate  region  is  wider  than  the  electron  inner  region  by  a  factor  of 
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1/4  1/4 

d./d„  ”  (T  /T  )  (m  /m  )  /2.  For  comparable  T. ,  and  T  ,  this 

1  e  li  el  i  e  il  el 

quantity  is  roughly  4.  The  intermediate  region  is  in  turn  thinner  than  the 
sheet  thickness  6  by  a  factor  of  /a^/26  <<  1.  In  this  paper,  we  treat  the 
three-region  matching  problem  (see  Fig.  2);  in  region  I,  the  total  current 
density  is  J1T  *  (J*"  +  J^)  ;  in  region  II,  J1T  =  (J*"  +  i  in  region 

III,  the  current  density  is  )  . 

B.  Linear  Dispersion  Relation 

The  linear  eigenvalue  equation  is  obtained  by  substituting  Eqs.  (22), 
(23),  (27)  and  (28)  into  Eq .  (4)  according  to  the  three-region  scheme 

described  in  the  preceding  section.  The  equation  is  then  solved  subject  to 
the  following  conditions.  The  solution  ij>  must  be  such  that  its  first 
derivative  Oty/dz)  vanishes  asymptotically  (|z|  +  ®)  and  that  the 
logarithmic  derivative  is  continuous  at  |z|  =  dg  and  |z|  =  d^.  For  the 

inner  region  (I),  the  eigenvalue  equation  is  found  to  be 


.  A2;(i) 

dz2 


(29) 


where  ^(x,z,t)  *  4>(z)  exp( ikx-iwt)  and 

Tii  ,1  6? 


A  ^  -  2ii  +  t^  mg-ng-)-  °i^ 

1  a2  i«  iH 

nr*  1 


)] 


+  r(i~+  MK 


kv 


ell 


en 


with  secli  (z/6)  ■  l  for  |z| 


<  dg  <<  6.  Here,  k  i  k6  and 


z/6. 


(30) 


By 


setting  a. 


a  *  0,  Eqs.  (29)  and  (30)  reduce  to  the  inner  equation  of 
6 


Ref.  8.  In  general,  A  is  independent  of  z  and  the  solution  is 


*ln(z)  *  C  cosh(Az) 


(31) 


where  C  is  an  integration  constant  and  the  symmetric  solution  has  been 
chosen. 

In  the  intermediate  region  II,  the  eigenvalue  equation  takes  on  the 

form 


•i 


i,  1  lkvi/J 
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[(kH2(kH  -  ‘e“(idHslnh'2<7)! 


el  e| 

A 

x  sech2(z)\|». 


el 


(32) 


It  is  easy  to  see  that  the  inner  ion  orbit  contribution  dominates  the  small 
gyroradlus  electron  contribution  in  this  region. 


In  the  exterior  region  |z|  >  d^  the  orbit  of  both  species  are  mainly 
of  the  small  gyroradius  drifting  type  and  the  eigenvalue  equation  takes  on 
the  form 


^  ■  [t2  -  2  ..=<.*<!>];  ♦ 

dz2  1  ill  il  ill 
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^kv 


el 


'•kv 


e| 


kv 


(33) 


el 


Asymptotically,  where  the  VB  contribution  vanishes  due  to  near  uniformity 
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of  the  equilibrium  magnetic  field,  the  solution  of  the  above  equation  is 


i|i  (z)  -  Dfl  +  —  tanhlzl  )expf-k|z| ),  (34) 

°°  k 

where  the  even  solution  is  chosen  and  D  is  an  integration  constant.  As 
required,  the  first  derivative  is  proportional  to  exp(-k|z|)  and  vanishes 
as  |  z  I  -*■  <». 

If  we  set  ”  a£  ■  0  in  Eq.  (30)  and  match  the  inner  and  asymptotic 
solutions  [Eqs.  (31)  and  (34)]  at  z  -  dg,  then  the  isotropic  results  of 
Ref.  8  are  recovered.  Equating  the  logarithmic  derivatives  of  Eqs.  (31) 
and  (34)  at  z  «  de,  we  obtain  the  linear  dispersion  relation 

A  tanh(Ade)  ■  (1  +  k  *tanh  dg)[k  ^sech^de  -  k  *(l  +  k  1tanh  de)], 

(35) 


where  A  is  given  by  Eq.  (30)  and  dg  =  d^/6 .  This  is  equivalent  to  the 
series  representation  in  terms  of  generalized  Legendre  functions  given  in 
Ref.  8.  Neglecting  terms  of  order  df  <<  1,  the  right  hand  side  can  be 
simplified  and  we  have  approximately 
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A  tanh  (Ad  )  -  —  -  k.  (36) 

e  k 

Using  the  small-argument  expansion  of  Z  function  in  Eq.  (32),  we  find 
(Ade)  «  1  for  the  low-frequency  perturbations.  Then,  keeping  only  the 
leading  terms  in  Eq.  (36),  we  find  that  the  Instability  Is  primarily  due  to 
the  resonant  electrons  and 

s2--£)3/2-^(>+r-)£-r).  07) 

Q 

which  is  identical  to  the  isotropic  result0.  It  shows  that  (y/kv  )  «  1 
2/2  _ 

scales  as  (ag/6)  and  that  y  >  0  (unstable)  for  k  <  1  and  y  <  0  (stable) 
for  k  >  1.  By  setting  y  m  0  and  neglecting  terms  of  order  d  ,  we  find 
the  marginal  condition 


0  ■*  —  -  k  , 
k 

o 

and  the  marginal  wavenumber  kQ  »  1,  recovering  the  previous  result  *  . 

For  the  anisotropic  case,  the  approximate  dispersion  relation  is  still 
given  by  Eq.  (36)  with  *  0  in  A  [Eq.  (30)].  As  we  can  see  from  Eq. 
(30),  the  electron  anisotropy  effects  are  large  because  of  the  small 
gyroradius  (a^/a^  «  m^/n^,  6  /ag  >>  1).  As  a  result,  it  is  likely  that 
nonlinear  effects  become  important  for  electrons  on  a  time  scale  comparable 
to  &>c^  This  makes  the  present  analysis  unsuitable  except  for  very 
small  degrees  of  electron  anisotropy,  viz.,  lTei/^e»  “  «  ae/5*  Thus, 

in  the  remainder  of  this  paper,  we  will  primarily  consider  the  case  with 
anisotropic  ions  and  isotropic  electrons.  Then,  keeping  only  the  leading 
terms,  we  have 

(Al.)2  •'  (1  +  ri)'1  +  2a,  (1  +  ^if l(f-)  ^  ,  (38) 

e  el  ell  ei  e  i 

where  -  (1  ~  T^/T^  ).  In  ®arginsl  3tate,  y  m  0  and 

(Ade  )2  -  20i(l  +!ii)'l(i-)  . 

el  e  i 

For  systems  in  which  (ae/6)  and  (m  /s^)  are  comparable  such  as  the 


09) 


_  2  _ 

magnetotail,  (Adg)  «  1.  Then,  the  marginal  wavenumber  kQ  is  given  by 

k  el  i  e 

o 

Using  Eqs.  (38)  and  (39)  in  Eq.  (36)  and  keeping  only  the  leading  terms,  we 
find 

e*-  - -j- £)3/2t>  -  £)[<i  - 1)  -  a 


!  II  /2-it 


ei 


el 


(40) 


This  is  the  approximate  anisotropic  dispersion  relation  obtained  by 
matching  the  inner  and  asymptotic  solutions  at  z  *  dg.  From  Eq.  (40),  we 
find  that  the  normalized  growth  rate  has  a  maximum  at 


and 


(r^trK1  +  ri)(>  +tc)- 

e  el 


(41) 


Equation  (36)  is  obtained  to  show  the  basic  properties  of  the 
isotropic  and  anisotropic  dispersion  relations.  One  noteworthy  property  of 
this  equation  is  that  A  is  either  purely  real  or  purely  imaginary.  Using 
small-argument  expansion  of  Z  and  W  functions,  it  is  easy  to  see  that  u  is 
nearly  purely  Imaginary  for  the  low-frequency  perturbations  being 
considered,  as  we  assumed  in  the  preceding  discussion  of  approximate 
solutions.  In  the  remainder  of  the  paper,  we  specialize  to  the  case  of 
imaginary  frequency  with  <«  =  iy  so  that  y  >  0  corresponds  to 
instability.  The  transcendental  equation  (36)  has  also  been  solved 
numerically  to  obtain  the  dispersion  relation  in  detail.  Figures  3  and  4 
show  the  growth  rate  y  normalized  to  the  asymptotic  ion  cyclotron 
frequency  versus  the  normalized  k  =  k6.  The  dispersion  curves  are 
shown  for  several  values  of  T^/T^  an<*  two  values  of  (aj/6),  0.03  and 
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4  Normalized  growth  rate  (y/u>cl)  versus  k6  for  a^/6  “  0.05  using 
Eq.  (37).  The  value  of  T^/Tj  is  (a)  0.9,  (b)  1.0,  (c)  1.1  and 
(d)  1.15.  (Tel/Te|  -  1). 


0.05.  In  particular,  the  value  a^/6  »  0.03  roughly  corresponds  to  the 
quiescent  magnetotail  parameters.  In  addition,  we  have  used  Tij/Tel  “  2. 
In  this  case,  ae/5  =  5  x  10  \ 

O 

The  curve  b  in  each  figure  corresponds  to  the  isotropic  result  with 

the  marginal  wavenumber  k5  “  1  independently  of  (a. /6).  Consistent  with 

13  1 

the  previous  conclusions  ,  we  find  that  temperature  anisotropy 

with  T^  /T  >  1  is  destabilizing  while  anisotropy  with  T^/T^  <  1  is 

stabilizing.  Here,  k_  is  parallel  to  so  that  the  case  with  the  greater 

temperature  perpendicular  to  k_  is  more  unstable.  These  figures  also  show 

that  Eq.  (39)  describes  the  stability  boundary  accurately.  For  example, 

for  T^/T  =  and  a^/6  “  0.03,  we  have  kQ  =  2.67.  For  a^/6  =  0.05,  we 

have  k  *  1.68,  in  agreement  with  these  figures.  In  addition,  Eq.  (40) 
o 

is  a  good  approximation  for  all  the  curves  shown  in  these  figures.  In 

comparing  Figs.  3  and  4,  note  that  the  two  values  of  a^/6  correspond  to 

different  values  of  6  and  ou  .  so  that  v  and  k  are  normalized  to  different 

cl 

scales  [see  Eqs.  (14),  (17)  and  (18)].  The  curves  are  not  completed  for 

k6  =  0  because  the  theory  breaks  down  as  k  +  0.  Comparing  curves  a,  c  and 

d  with  curve  b  in  each  of  Figs.  3  and  4,  we  see  that  the  effect  of 

anisotropy  is  substantial.  This  is  expected  because  the  anisotropy  term, 

the  third  term  in  Eq.  (8),  is  greater  than  the  isotropic  term  by  a  factor 

1  2 

of  (kv,  /a))  »  1.  In  this  regard,  we  point  out  that  Laval  and  Pellat 
J  ^ 

used  an  energy  method  to  show  that  the  mode  treated  here  is  completely 
stabilized  for 


(42) 


For  the  parameter  used  in  Fig.  3,  the  right  hand  side  is  approximately 
0.9995.  Using  this  value  of  electron  temperature  anisotropy  in  Eq.  (37), 
we  find  that  the  mode  is  in  fact  completely  stabilized.  With  electron 
anisotropy,  the  square  brackets  in  Eq.  (40)  would  contain  another  term 
identical  to  the  second  term  except  for  the  replacement  e  ~  i  and 
me/m1  1.  [See  the  discussion  preceding  Eq,  (38)].  The  reason  for  the 
extremely  sensitive  dependence  on  the  electron  anisotropy  described  above 
is  that  the  electron  term  without  the  small  mass  ratio  enhances  the 
anisotropy  effects  discussed  above.  As  a  general  remark,  we  note 
that  Y/«ci  18  typically  of  the  order  of  10”^.  In  addition,  y/kv^  is  also 
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of  the  order  of  10-4  for  the  unstable  parameter  regimes  for  both  species. 
Thus,  the  low  frequency  approximation  used  to  derive  Eq.  (8)  is  justified 
posteriori. 

So  far,  we  have  examined  anisotropic  tearing-mode  properties  using  the 

two-region  matching  method  following  a  number  of  previous  works.  However, 

examination  of  Eq.  (32)  shows  that,  in  the  intermediate  region,  it  is  the 

ion  orbits  that  dominate.  Furthermore,  because  d^  >>  dg,  the  effect  of 

these  ions  is  expected  to  be  large.  In  order  to  study  the  anisotropic 

properties  including  the  intermediate  region,  we  have  numerically 

integrated  Eqs.  (32)  and  (33)  in  regions  II  and  III.  The  physically 

acceptable  solution  must  have  the  asymptotic  behavior  given  by  Eq.  (34)  and 

the  logarithmic  derivative  must  be  continuous  at  z  =  d^  and  z  =  dg.  At  z  = 

de,  the  derivative  is  matched  to  that  of  the  analytic  solution,  Eq.  (31). 

The  matching  condition  then  gives  the  linear  anisotropic  dispersion 

relation.  The  results  are  plotted  in  Figs.  5  and  6  for  several  values 

of  T.  /T  and  for  two  values  of  a. /6,  0.03  and  0.05.  Here,  the  electrons 

ii  i  II  i 

are  isotropic. 

Comparing  Figs.  (3)  with  (5)  and  (4)  with  (6),  we  see  that  inclusion 
of  the  ion  orbits  in  the  intermediate  region  modifies  the  growth  rate  and 
the  stability  boundary  significantly.  In  general,  the  anisotropy  effects 

are  further  enhanced  by  the  inclusion  of  the  intermediate  region  (region 

II).  That  is,  for  T^/T^  ^  1*  th e  mode  is  more  strongly  stabilized  with 
the  intermediate  region  than  without  it.  For  T^  /T^  >  1,  the 

instability  is  more  strongly  enhanced  with  the  intermediate  region  than 
without  it.  However,  the  isotropic  dispersion  curves  obtained  using  the 
"three-region”  matching  method  are  nearly  identical  to  the  corresponding 
ones  obtained  using  the  two-region  approach.  The  absolute  value  of 
of  curve  b  in  Fig.  3  is  greater  than  that  of  curve  b  in  Fig.  5  by 
approximately  3%  to  4%  (a^/6  -  0.03).  The  absolute  value  of  (y/u)^)  of 
curve  b  in  Fig.  4  is  also  greater  than  that  of  curve  b  in  Fig.  6  by  similar 
amounts  (a^/6  ■  0.05).  This  agreement  is  nontrivial  since  Eq.  (32)  shows 
that  the  dominant  ion  orbits  in  the  intermediate  region  modifies  the 

j 

eigenvalue  equation  significantly  and  indicates  that  the  isotropic 
dispersion  relation  is  indeed  determined  primarily  by  the  resonant 
electrons  in  the  inner  region  |z|  <  dfi.  The  above  behavior  can  be 
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Normalized  growth  rate  (y/<ucj)  versus  kS  for  a ^/6  *  0.03 
including  the  intermediate  ion  region.  The  value  of  /T 
(a)  0.9,  (b)  1.0,  (c)  1.05,  (d)  1.1  and  (e)  1.15.  (T  /T 


understood  from  Eqs.  (6)  or  (8)  by  noting  that  ikqSy  is  the  time-integrated 
Lorentz  force  in  the  x-direction  causing  the  particles  to  bunch.  Since  the 
ion  orbits  are  larger  than  the  electron  orbits  by  approximately  d^/d^  4, 
the  accumulated  influence  is  also  increased.  It  is  evident  that  the 

Lorentz  force  term  is  zero  in  the  isotropic  case.  It  is  of  interest  to 

note  that  the  force  responsible  for  the  anisotropic  effects  is  similar  to 

20 

that  in  the  mirror  instability  .  However,  unlike  the  mirror  instability, 
there  is  no  threshold  value  of  T^/T^  that  is  required  for  the  onset  of 
the  anisotropic  effects. 

Figures  5  and  6  show  that  for  T  /T^  >  l  the  maximum  growth  rate  is 

enhanced  by  one  order  of  magnitude  or  more  from  the  isotropic  case  and  that 

the  marginal  wavenumber  is  increased.  For  T^/T  ^  the  instability  is 
essentially  stabilized.  Note  that,  as  before,  (y/wcl)  is  typically  of  the 
order  of  10_l!t  so  that  the  low  frequency  approximation  is  justified.  The 
results  are  shown  for  /T^  up  to  1.15.  For  higher  >  1.3) 

degrees  of  anisotropy,  the  increasing  value  of  (y/w^)  would  render  the  low 
frequency  and  constant-^  approximations  invalid.  The  dashed  line  (e)  in 
Fig.  6  shows  the  dispersion  curve  for  T^  /T^  =  1.1  with  the  VB 
contribution  neglected.  We  see  that  the  qualitative  behavior  is  not 
significantly  changed  and  that  the  VB  drift  has  the  opposite,  albeit  small, 

effect  from  the  axis-crossing  orbits.  This  is  easy  to  understand  siftfc  the 

guiding  center  of  a  drifting  particle  is  opposite  to  the  mean  drift 
velocity  of  the  plasma  layer  [Eq.  (25)]. 

Figure  7  shows  the  maximum  growth  rate  (y/j  .)  as  a  function  of 

ci  max 

temperature  anisotropy  (T,  /T.  ).  In  Fig.  8,  we  have  plotted  the  eigen- 

A  a  X  J-  II 

function  4>(z)  for  |z|  >  dg  for  two  values  of  T^/T_^ .  Both  curves 
correspond  to  the  respective  maximum  growth  rates  (Y/<*>ci)max»  The  matching 
surfaces  at  Izl  *  dfi  and  izi  ■  d^  are  also  shown.  We  see  that  the 
constant-41  approximation  is  reasonable  for  the  values  of  T^/T^  used 
here.  However,  for  T^/T^  >  1.3,  the  variation  in  41  within  the  ion  layer 
becomes  substantially  greater. 
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Normalized  growth  rate  (Y/wci)  versus  k5  for  a^/6  *  0.05 
including  the  intermediate  ion  region.  The  value  of  T^/T^  is 
(a)  0.9,  (b)  1.0,  (c)  1.05,  (d)  1.1,  (e)  1.1  with  VB  neglected 
and  (f)  1.15.  (Tei/Te|  -  1). 


IV.  SUMMARY 


In  this  paper,  we  have  investigated  the  collisionless  tearing-mode 

properties  of  an  anisotropic  (T^  *  T ( )  neutral  sheet.  Both  ions  and 

electrons  are  described  by  Vlasov  distribution  functions.  Using  simplified 

particle  orbits  and  constant-41  approximation,  the  perturbed  distribution 

function  is  evaluated  for  low  frequency  perturbations  (y/<»>c1  «  1)  with 

k  n  B  .  The  first-order  current  densities  are  explicitly  found  [Eqs. 
—  — o 

(22),  (23),  (27)  and  (28)]  and  eigenvalue  equation  is  obtained  for  the 

three  regions  [Fig-  2  and  Eqs.  (29),  (32)  and  (33)].  The  equation  is  solved 

using  both  analytic  approximations  and  numerical  methods  to  obtain  the 

linear  dispersion  relation  (Figs.  3,  4,  5  and  6)  and  the  eigenmode 

structure  (Fig.  8).  First,  by  neglecting  the  ion  intermediate  region 

(region  II),  the  conventional  two-region  matching  method  Is  used  to  find 

the  approximate  anisotropic  dispersion  relation  [Eq.  (40)  and  Figs.  3  and 

4]  as  well  as  other  dispersion  properties  [Eqs.  (39)  and  (41)].  Then,  the 

full  eigenvalue  equation  is  numerically  integrated  in  regions  II  and  III. 

The  three-region  matching  condition  then  gives  rise  to  the  anisotropic 

dispersion  relation  illustrated  in  Figs.  5  and  6  for  a  number  of  parameter 

values.  It  is  shown  that  temperature  anisotropy  with  T,  /T.  >  1 

i±  m 

enhances  the  growth  rate  by  as  much  as  an  order  of  magnitude  or  more  while 

anisotropy  with  T  /T  ^  *  strongly  stabilizes  the  mode.  This  is 

•L  « 

consistent  with  previous  results  based  on  energy  principle  and  marginal 
13 

stability  calculations.  It  is  also  found  that  the  conventional  approach 
of  matching  the  inner  and  outer  solutions  at  the  electron  inner  region  (|z! 
<  dg)  is  not  adequate  in  the  anisotropic  case.  An  intermediate  region 
(de  <  I z |  <  d^)  is  identified  in  which  the  axis-crossing  ion  orbits  are 

dominant.  It  is  the  ions  in  this  region  that  account  for  the  substantial 
differences  as  exhibited  by  Figs.  3,  4,  5  and  6.  The  use  of  simplified 
orbits  also  allows  evaluation  of  the  effects  of  different  orbits 
explicitly.  In  particular,  it  is  found  that  the  inner  orbits  and  outer 
>rblts  (the  VB  drift  orbits)  have  the  opposite  effects  on  the  growth 
rate:  where  the  inner  orbits  are  destabilizing,  the  outer  orbits  are 

stabilizing  and  vice  versa  In  the  anisotropic  case. 


Eigenfunction  ij»(z)  versus  z/6  corresponding  to  (v<*>cl)max  fo 


(a)  T±  /T  *  1.05  and  (b)  T^/T^  -  1.1.  -  0.03.  The 

matching  boundaries  are  indicated:  |zl  3  dg  and  |zl  »  d^. 

(T.i/T..  '  '>• 


A  physical  system  for  which  the  present  analysis  may  be  applicable  is 
the  earth's  magnetotail.  The  previous  works  on  the  tearing-mode  stability 
of  such  a  system  have  studied  isotropic  neutral  sheets  (see,  for  example. 
Refs.  8,  10  and  11).  In  light  of  the  fact  that  the  magnetotail  is  highly 
collisionless,  it  is  reasonable  to  expect  the  temperature  distribution  to 
be  generally  anisotropic.  Our  present  results  suggest  that  the  linear 
tearing-mode  properties  of  the  magnetotail  and  similar  physical  systems  are 
dominated  by  the  anisotropic  tearing-mode.  In  fact,  the  It  n  mode  can 
grow  significantly  faster  than  previously  predicted  if  temperature 
anisotropy  is  present  (c^  >  a^S).  Thus,  the  particle  temperature 
distribution  is  a  critically  important  quantity  for  understanding  the 
linear  tearing-mode  stability  properties. 

The  present  analysis  utilized  simple  approximate  orbits  in  evaluating 

Q 

the  orbit  integrals.  In  the  isotropic  case,  the  two-region  approach 
yielded  results  in  agreement  with  numerical  results,^  indicating  that  the 
stability  properties  are  not  critically  dependent  on  the  precise  orbits. 
In  the  anisotropic  case,  however,  the  large  ion  orbits  are  important  so 
that  a  more  accurate  calculation  of  the  orbit  integrals  is  desirable.  In 
the  present  analysis,  the  orbits  that  are  neglected  are  mainly  those  of 
particles  with  energy  substantially  greater  than  the  typical  velocities, 
constituting  a  relatively  small  fraction  of  the  total  particles.  In 
conjunction  with  the  fact  that  the  low  frequency  and  constant-^  approxima¬ 
tions  are  well  satisfied  by  the  results,  we  expect  the  present  results  to 
be  a  good  first  approximation  for  understanding  the  essential  physics  of 
the  collisionless  anisotropic  tearing-mode  properties. 
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